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Abstract. The matrix-valued Wcyl-Titchmarsh functions M(A) of vector-valued 
Sturm-Liouville operators on the unit interval with the Dirichlet boundary conditions 
are considered. The collection of the eigenvalues (i.e., poles of M(A)) and the residues 
j^-V) ■ of A/(A) is called the spectral data of the operator. The complete characterization 

, of spectral data (or, equivalently, N x N Weyl-Titchmarsh functions) corresponding 

q^' to N x N self-adjoint square-integrable matrix-valued potentials is given, if all N 

, eigenvalues of the averaged potential are distinct. 

m ; 

JT" 1 ! 1. Introduction 

Ph . 

^ \ We start with a short description of known results in the inverse spectral theory 

for scalar Strum-Liouville operators on a finite interval. We recall only some important 
steps mostly focusing on the characterization problem, i.e., the complete description 
of spectral data that correspond to some fixed class of potentials. More information 
about different approaches to inverse spectral problems can be found in the monographs 
(N '■ [Mar86] . |Lev87j . |PT87j . jFYOlj . survey [Ges07] and reference s there in. 

^ . The inverse spectral theory goes back to the seminal paper p3o46| (see also |Le49j ). 

■^j- ! Borg showed that spectra of two Sturm-Liouville problems —y" + q(x)y = \y, x G [0, 1], 

with the same boundary conditions at 1 but different boundary conditions at 0, deter- 
mine the potential q(x) and the boundary conditions uniquely. Later on, Marchenko 
[Mar 50] proved that the so-called spectral function p(A) (or, equivalently, the Weyl- 
Titchmarsh function m(A)) determines the potential uniquely. Note that the spectral 
function is piecewise-linear outside the spectrum {A™}^ and its jump at A n is equal 
to the so-called normalizing constant [a n (q)]~ l given by (11.31) . At the same time, a 
^ different approach to this problem was developed by Krein [Kr51j . |Kr53j . [Kr54] . 

An important result was obtained by Gel'fand and Levitan |GL51j . They gave 
an effective method to reconstruct the potential q from its spectral function. More 
precisely, they derived an integral equation and expressed q(x) explicitly in terms of 
the solution of this equation. At that time, there was some gap between necessary and 
sufficient conditions for the spectral functions corresponding to fixed classes of q(x). 
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Some characterization of spectral data for q such that q^ m ' G >C 1 (0, 1) was derived 
by Levitan and Gasymov [LG64J for all m = 0,1,2,... Also, they gave the solution 
of the characterization problem in the case q" G £ 2 (0, 1). Marchenko and Ostrovski 
|MU75j obtained a sharpening of this result. Namely, for all m — 0, 1, 2, .. they gave 
the complete solution of the inverse problem in terms of two spectra, if q^ G £ 2 (0, 1). 

Trubowitz and co-authors (Isaacson [IT88] . McKean |IMT84j . Dahlberg [DT84j . 
Poschel [PT87J) suggested another approach. It is based on the analytic properties of 
the mapping {potentials} t— > {spectral data} and the explicit transforms corresponding 
to the change of only a finite number of spectral parameters (A n (g), v n (<l))t=v Their 
norming constants v n {q) differ slightly from the normalizing constants (ll.3p . but the 
characterizations are equivalent (see Appendix [B]) . Also, this approach was applied to 
other scalar inverse problems with purely discrete spectrum (singular Sturm-Liouville 
operator on [0, 1] [GR88] : perturbed harmonic oscillator |MT81| . [CKK04] . [CK07]). 

Thus, nowadays the inverse spectral theory for the scalar Sturm-Liouville opera- 
tors is well understood. By contrast, until recently only some particular results were 
known for vector-valued operators. 

In our paper we consider the inverse problem for the self-adjoint operators 

L^ = -< + V(x)V, V(o) = = o, Ve£ a ([o,i];C"), (1.1) 

where V = V* G £ 2 ([0, 1]; C NxN ) is a self-adjoint N x N matrix-valued potential. 
Denote by <f(x) = <p(x, A, V) and x( x ) = x( x -> \ V) the matrix-valued solutions of the 
equation — ip" + V(x)ip = Xip such that 

^(0) = x(l) = 0, y/(0) = = I N , 

here and below In denotes the identity NxN matrix. Note that 

x {x,A,V) = yp(l-x,A,1/ s ), where V\x) = V(l — x), x G [0, 1]. 

The matrix-valued Weyl- Titchmarsh function for this problem is given by 

M(A) = M(A, V) = [x'y-'KO, A, V) = [M(A)]*, A G C. (1.2) 

In the scalar case, the Weyl-Titchmarsh function m(A, q) is a meromorphic function 
having simple poles at Dirichlet eigenvalues X n (q) and 



res m(X,q) = -[a n (q)}' 
X=X n (q) 



\<p(x, \ n ,q)\ 2 dx 



;i.s) 



So, the sharp characterization of all scalar Weyl-Titchmarsh functions (or, equiva- 
lently, all spectral data (A n (g), 0£ n (q))^ 1 ) that correspond to potentials q G £ 2 (0, 1) is 
available due to [M075J or [PT87J (see also Appendix [B]). Namely, the necessary and 
sufficient conditions are 

Ai < A 2 < A 3 < (A n - Tc 2 n 2 - q )n=i ^ ^ 2 f° r some 9o ^ M 

and {nn ■ (2n 2 n 2 a n (q) - 1))+^ G £ 2 . ' ' " ' ' 

In the vector-valued case, it is known that the Weyl-Titchmarsh function deter- 
mines V uniquely (see [Mal05j or [Yur06j ) . Some other miscellaneous results concern- 
ing vector- valued Schrodinger operators were obtained in |Car02] . |CK06aj . [ChSh97j . 
[CHGLQ0], [.TL98aj . [,TL98bj . [SFHi] . [ShnT] . Nevertheless, to the best of our knowledge, 
no solutions of the characterization problems have been available until recently. 
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Following [CK06b] . we denote by Ai < A 2 < .. < X a < ... the eigenvalues of 
L and by k a = dimE a G [1, iV] their multiplicities, where E Q C £ 2 ([0, 1]; C^) is 
the eigenspace corresponding to the eigenvalue A a . Then (see details in |CK 06b|). 
the Weyl-Titchmarsh function M(A) is meromorphic outside the Dirichlet spectrum 
a(V) = {X a (V)} a>l and 

res M(A) = -B a = -p^g^Pa, 

where 

p a : C N - £ a = Kerp(l, X a ,V) = {heC N : ^ = </?(•, A a , V)h e E a } 
is the orthogonal projector and 

g a =Pa I [(p*(p](x, X a , V)dx p* a = g* a > 



[(p*(p](x, X a , V)dx 
Uo 

is the self-adjoint operator (or the normalizing matrix) acting in £ a . We also use the 
notation P a = p* a p a : — > £ a C C^. Note that for all hi, hi G £ a one has 

( < 0a;/n,'0a;/ l 2)£ 2 ([o,i];c^) = / K[(p*(p] (x, X a , V)hi dx = {hi, g a h 2 )e a ■ 

Jo 

We call (A a , P a , ga)a=i the spectral data of the operator L. If k a = 1, then g a acts in 
the one-dimensional space £ a , so we consider it as a positive real number (and call it, 
as in the scalar case, the normalizing constant). The spectral data determine (e.g., see 
Proposition 12.61) the function M(A), and so the potential V(x), uniquely. The main 
result of our paper is the following solution of the characterization problem. 

Let e\,e^, be the standard coordinate basis and P® — (')^)e° be the coor- 
dinate projectors in C . We denote the Euclidian norm of vectors h G and the 
operator norm of matrices A G <C NxN by \h\ and \A\, respectively. 

Theorem 1.1 (Characterization of spectral data). For all v® < u° < ■■ < the 

mapping V i— > (A a , P a , g a )ti=i ^ s a bijection between the space of potentials 

V = V* G C 2 ([0, 1]; C NxN ) suchthat [ V(x)dx = diag{u?, v%, .., v%} (1.5) 

Jo 

and the class of spectral data satisfying the following conditions (A)-(C): 

(A) The spectrum is asymptotically simple, i.e., there exist a ^ ; n° ^ 1 such that 

kl + kt + .. + Ko =JV(n°-l) and k° a = 1 for all a ^ a°+l. 

It allows us to define the double-indexing (n,j), n^n°, j — 1,2, ..,N, instead of 'a>a°. 
Namely, we set X n j = X a o + N{ n ~n Q )+j, Pn,j — Pa°+N(n~n < >)+j and so on for n^n® . 

(B) The following hold true for all j = 1, 2, .., N: 

(A nj -7r 2 n 2 -^)+r n o G t 2 , (tto • {2it 2 n 2 g n j-l))+™ n « G £ 2 , 

(|P ni ,-P°|)+- G f and (vm • |£f =1 P nj -/;v|& G £ 2 . ' 

(C) 27ie collection (X a ; P a )„^ satisfies the following property: 

Let £ : C — > C N be an entire vector-valued function. If P a £{X a ) = for all 
a^l, f(A) =0(e |ImVX| ) as |A| -> oo and £ G £ 2 (M+), taen £(A) = 0. 
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Remark 1.2. Let V — V* G £ 2 ([0, 1]; C NxN ). Applying some unitary transform in C N , 
one may always assume that V(x)dx = diag{t> J 1 , t> °, --^n}' v i ^ v 2 ^ •• ^ v n- Our 
assumption ( li.5j) states that all the v® are distinct. It simplifies the analysis, since 
otherwise infinitely many eigenvalues X a can be multiple. In particular, in the general 
case, one has to introduce some other parameters instead of(P n j, g n ,j)- 

We give also a simple reformulation of the algebraic restriction (C) (note that it 
doesn't depend on the shift of the spectrum). 

Proposition 1.3 (reformulation of (C)). Let X a > for all a ^ 1 and P a = h a h* a , 
where h a = (ha^ ; .. ; h^a a ^) consists of k a orthonormal vectors ha G C N . Then the 
condition (C) is equivalent to the following: 

Vector-valued functions e ±% '^" t h < £ > , j = 1, .., k a , a ^ 1, together with the 
constant vectors e?, .., e° N span £ 2 ([— 1, 1] ; C N ). 

Remark 1.4. In the scalar case, (C) always holds true due to the well known result of 
Paley and Wiener (e.g., see |Le40j p.J^l). In the vector-valued case, this condition is 
not trivial. Some discussion of ( C) is given in AppendixlAl (see Propositions \A.3l \A.4\ ). 
Note that, if P n ,j = Pj? for all n ^ m + 1 and j = 1, 2, .., N, then one can reformulate 
(C) as the condition det T ^ for some NmxNm matrix T (see Proposition \A.5\) . 

As usual, Theorem 11.11 consists of several different parts: 

(i) Uniqueness Theorem (spectral data determine the potential uniquely); 

(ii) Direct Problem (spectral data constructed by a given potential satisfy (A)-(C)); 

(iii) Surjection (any data satisfying (A)-(C) are spectral data of some potential). 

We do not discuss the uniqueness theorem (i) in our paper and refer to [Mal05j . 
|Yur06j (or [CK06b] ) for this fact. The direct problem (ii) is considered in Sect.HJ Note 
that the spectrum is asymptotically simple due to our assumption t>° < < .. < 
(see also Remark 11.2ft . As in the scalar case, the Fourier coefficients of V appear as 
leading terms in the asymptotics of the spectral data (Propositions [27T1 and I2~5l) . We 
also give the explicit expression for M(A) in terms of the spectral data in Sect. 12.41 

The main part of our paper (Sect. [3]) is devoted to the surjection (iii). The 
general strategy of the proof is described in detail in Sect. 13 . 1 1 Here we give only a 
short sketch of our arguments. We start with some admissible data (A*, Pt, g a ) a ^i 
satisfying (A)-(C). Using the well known characterization ( 11.41) for the scalar case, we 
construct some special diagonal potential V° such that ^(U ) = {A^} a ^i. 

In Sect. l3.2H3.4l we introduce some essential modification of the spectral data in or- 
der (a) to control the splitting of multiple eigenvalues and (b) to join together all asymp- 
totics in ( 11. 6p . We prove that the mapping $ : {potentials} i— >■ {modified spectral data} 
is real-analyticQ near V°. The main purpose of involving analyticity arguments here is 
the well known equivalence of the analyticity and the weak-analyticitjo for mappings 
between complex Hilbert spaces. Thus, we immediately derive the smoothness of the 
whole mapping $ from the smoothness of its components. 



1 The mapping F : U — > H ^ between real Hilbert spaces U C and is real- analytic iff it 
has continuation Fc ■ Uc — * into some complex neighborhood U Ct/cC that is differ entiable 
as the mapping between the complexifications \ of the real spaces H^K 

2 In Hilbert spaces, the weak-analyticity is equivalent to the analyticity of particular coordinates 
and the local boundedness, see nice Appendix A in [PT87J or the monograph |Di99j for details. 
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In Sect. 13.51 13.61 we use the Fredholm Alternative in order to show that $ is a 
local isomorphism near V (i.e., <iyo$ is invertible). Thus, all additional spectral data 
sufficiently close to (P a (V ), (? a (V r<> )) Q ^i can be obtained from potentials having the 
same spectrum {A*} a ^i as V°. In particular, if a* is large enough, then there exists 
V such that a(V) = {K)^i and (P Q (V) , g Q {V)) = {Pt,g ] a ) for all a > a'. 

We complete the proof in Sect. 13.71 using the explicit isospectral transforms con- 
structed in our recent paper [CKQ6bj . As usual in Trubowitz's approach, we need to 
change only some finite number a* of additional spectral data {P a ,g a )- Note that the 
condition (C) and the restrictions introduced in [CK06b] in terms of "forbidden" sub- 
spaces are equivalent (see Proposition IA.4I) . Thus, one can change any finite number 
of projectors P a in an arbitrary way that doesn't violate (C) (see details in Sect. 13 . 71) . 

Note that we do not present any explicit reconstruction procedure for the potential, 
if there are infinitely many perturbed spectral data. The natural idea is to use some 
passage to the limit changing the residues BaiV ) i— > B^, a = 1,2,.., of the Weyl- 
Titchmarsh function step by step. Each step is doable due to isospectral transforms 
constructed in [CK06bJ but we do not prove the convergence of this procedure. 

We finish the introduction with several remarks concerning some possible further 
developments of our approach to this inverse problem. 

Remark 1.5. The isospectral transforms constructed in |CK06b| generalize the scalar 
isospectral flows (see [PT87] ) and some specific class of isospectral transforms given in 
[JL98aj . Nevertheless, to the best of our knowledge, no analogues of the explicit flows 
changing the eigenvalues (see |PT87j ) are known in the vector- valued case. We think 
that such a construction would simplify the inverse theory a lot. 

Remark 1.6. One may be interested in the characterization for other parameters, e.g. 
the spectra of several boundary problems (similarly to the original paper [Bo46j). 
Almost nothing is known here. Yurko [Yur06j proved that A^ 2 + l spectra determine 
the potential uniquely. On the other hand, the naive count says that this inverse prob- 
lem is overdetermined. Note that, in the spirit of Appendix [Bj this question can be 
considered as a parametrization problem for some class of matrix- valued functions. 

Remark 1.7. Consider the Schrodinger operator Hy = —y" + Vy on R with a. NxN 
potential V = V* such that |:r|)|V(:r)|Gfe < +oo (e.g., see |Q185j ). It has a finite 

number of eigenvalues Ai < .. < A m < with the multiplicities k a = dimEc, where E Q 
is the eigenspace corresponding to X a . In order to solve the inverse scattering prob- 
lem completely, one needs to characterize the residues of the transmission coefficient 
at X a . Unfortunately, we do not know any results in this direction. For the scattering 
problem on the half-line a characterization was given in [AM63] but it involves implicit 
conditions for spectral data (much more complicated than our condition (C)). 

Remark 1.8. In the scalar case, the Dirichlet eigenvalues and the norming constants are 
canonically conjugate variables for the Korteweg-de Vries equation with periodic initial 
conditions (see [FM76J). Similarly, the (negative) eigenvalues and the corresponding 
normalizing constants of the (scalar) Schrodinger operator — y" + q(x)y on IR with a 
decreasing potential q(x) are canonically conjugate variables for the Korteweg-de Vries 
equation (see [ZF71] ) . The vector- valued case is more complicated (see [CD76j , [CD77j , 
[0185]). We hope that our results could be useful from this point of view. 
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2. Direct problem 
2.1. Asymptotics of the eigenvalues and the individual projectors. Denote by 

t/(°) = ! V(t)dt, V {cn) = [ V(t) cos 2-rmtdt and V {sn) = [ V(t) sin27rnt dt 
Jo Jo Jo 

the (matrix) Fourier coefficients of V. We start with some elementary asymptotics 

of the fundamental solutions ip(x, A, V) and A, V) = ip(l — x, A, V") for A close 

to ir 2 n 2 . It's well known that 

sinzx 1 f x / e \lmz\x\ 

ip(x, z 2 , V) = I N + — sin z(x-t) ■ V(t) sin ztdt + Ol , .„ . (2.1) 

z z Jo \ \ z \ J 

Here and below constants in O-type estimates depend on the potential. In this section 

we do not pay the attention to the nature of this dependence. Let 

2 2 2,.. .. , A* , r\l ^ 



z = ir n + a, a = 0(1), so z = nn H hO, 

2nn \n 6 

Then, 

sin zx I f x 

ld(x,z 2 ,V) = In H — / smirn(x— t) ■ V(t) sinirntdt + O 

nn 7i 2 n 2 J 

In particular, 

_iy» 



n 3 



tp{l,z 2 ,V) 



27c 2 n 2 



fiI N -V^ + V^ + o(^j 



(2.2) 



Proposition 2.1. Let V = V* e £ 2 ([0, 1]; C NxN ) satisfy V<® = diag{< <..,<} 
with < v® < •• < v N . Then, 

(i) there exists n° = rf(V) ^ ||V|| such that (a) there are exactly N(n°—1) eigenvalues 
counting with multiplicities in the interval (— 7r 2 (n° — l) 2 — 3||V||; 7r 2 (n° — l) 2 + 3||l / ||), 
(b) for each n ^ n° there are exactly N simple eigenvalues \ n> i, A„ )2 , .., A ni jv in the 
interval (7i 2 n 2 — 3||V||; 7r 2 n 2 + 3||V||), (c) there are no other eigenvalues; 

(ii) for each j = 1, 2, .., N the following asymptotics hold true as n — > oo: 

Xn . = k 2 u 2 + v° - tf£> + 0(S n (V)), where S n (V) = \V^\ 2 + -; 



n 



(Hi) if p n ,j = {^h n j)h nt j, where h n j G is such that \h n j\ = 1, (h n> j, e°) > ; then 
the asymptotics 

^4cn) ^4cn) ^4cn) ^~(cn) 



\v%- ^ u „ 1 - v u j+1 — v)j v N - v]j ' 

hold true for each j = 1,2, .., N as n — > oo. 
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Note that the condition n (V) ^ ||V|| guarantees that the mentioned intervals do not 
intersect each other. We need the following simple matrix version of Rouche's Theorem: 



Lemma 2.2. Let F,G : B(w,r) — > C be analytic matrix-valued functions such that 
|G(A)| • |F _1 (A)| < 1 for all A on the boundary of some disc B(w,r) C C. Then, 
the scalar functions det F and det(F + G) have the same number of zeros in B(w,r) 
counting with multiplicities. 

Proof. We check that Ac arg(det F) = Ac arg(det(F + G)), where A^arg/ denotes 
the increment of arg/ along the circumference C = {A : |A — w\ = r}. Note that, if 
A G C, then all eigenvalues of I + G(X)F^ 1 (X) have strictly positive real parts since 
|G(A)F- 1 (A)| < 1. Thus, the result follows from 

A c arg(det(F+G)) - A c arg(det F) = A c arg(det(/+GF" 1 )) = 

and the classical argument principle. □ 

Proof of Proposition \2.1[ (i) Firstly, we apply Lemma 12.21 to the function 

x(0,\,V)=<p(l,\,V*) = F(\)+G(\) 

in the discs 

{A : |A| <7rV+3||V||} with F(X) = ^^I N 

v A 

(see asymptotics (12.11) ) and 

{A: A = vrV + /i, \n\ <3\\V\\} with F(X) = ^44 ( (A-ttV)/ - V {0) ) 

(see asymptotics (I2.2p ). Thus, if n is sufficiently large, then there are exactly Nn and 
N eigenvalues (zeros of det x(0, •, V)), respectively, inside these discs counting with 
multiplicities. Secondly, let 

d = \ min i= i v . iAr _i(f° +1 - vf). 

If n is sufficiently large, then |T/( cra )| is small and one can apply Lemma [2.21 (with the 
same functions F as above) in the discs 

{A : A = vrV + v° + n, < d}, j = 1, 2, .., N. 

So, if n ^ n , then there are exactly one simple eigenvalue X n j = n 2 n 2 + fi n j inside 
each small disc B(n 2 n 2 +Vj,d) and there are no other eigenvalues. 

(ii) Recall that det ip(l, X n j, V) = 0. Therefore, due to (12.21) and the standard pertur- 
bation theory, the self-adjoint matrix Unjl^ — 

y(0) + y(cn) hag at leagt 

one eigenvalue r 

such that \t\ — 0(n _1 ). On the other hand, the eigenvalues of the matrix — 
are t s = v ° — viT^ + 0(\V( cn ^\ 2 ), s = 1, 2, .., N. Hence, for some s, 

Vn,j - v° s + u<f ) = 0(\V^\ 2 ) + Oin- 1 ). 

Due to (i), s — j. 



8 



DMITRY CHELKAK AND EVGENY KOROTYAEV 



(iii) Let j — 1 for the simplicity and d° k — v® — v k , k = 2, .., N. In view of (12. 2p and (ii) 



Men) JO _ Men) , Men) Men) 
u 21 a 2 u ll > u 22 ■■■ U 2N 



2ir 2 n 2 



o 



Sn(V) 

n 2 



\ tfen) Men) , _ Men), Men) I 

\ u Nl U 2N ••• a N u ll ' v NN / 

Recall that <p(l, A n> i, V)h Uj \ = 0. Thus, 

A n>1 , V)h nA , el) = gives (h nA , e° k ) = 0(\V^\ + 5 n (V)) for all k = 2, .., N, 
\h nA \ = 1 gives (h nA , e?) = 1 + 0{5 n {V)) 
and, using ((p(l, X n> i, V)h n> i, e° k ) = again, one obtains 

4? } + 4 ■ (Ki, e °) + 0(5 n (V)) = 0, k = 2, .., N. 
Note that (ii), (iii) are standard results for the perturbation of a simple eigenvalue. □ 

2.2. Asymptotics of the norming constants and the averaged projectors. 

Due to Proposition 12.11 all sufficiently large eigenvalues are simple. Therefore, for all 
sufficiently large n ^ n° and j = 1, 2, .., N we may introduce the factorization 

P n ,j = h n ,jh* nJ , B n>j = - res M(A) = h n>j g^h nJ = g~jP n>j , 
where g n j > 0, h n j G C N , \h n j\ = 1 and (/t n j,e°) > 0. Denote 

N 

B n = B n (V) = Y,B n ,j, n>n\ 
i=i 

We begin with some simple reformulations of the needed asymptotics. Note that 
Proposition 12.11 gives 

h nJ = e° + f for all j = 1, 2, .., N. (2.3) 
Here and below we write a n = b n + £^ iff 

(l a n — b n \)n= n o & ll = {(c n )n= n o ■ (n k c n )n= n o € £ 2 } . 
Lemma 2.3. The following asymptotics are equivalent: 

(ii) (/i nJ , /i„, fc ) = /or a// j ^ k, j, k = 1, 2, .., N. 
Proof. Introduce N x N matrices h n = ( h nj ± ; h Hj 2 ; ... ; h n ^ ). Then 

N N 

h n h n = ^ ] ^ l n,jh n j = ^ ] P n j 

and 

^n^n = (h* n ,j h n,k) j k=1 = {{K,k, K,j)) j k=l ■ 

The matrices h n h* n and are unitary equivalent (since h n h* n = u n (h n h n )u n , where 
the polar decomposition of h n ). Thus, the asymptotics h n h* n = I^ + if are 
equivalent to the asymptotics h* n h n = 1^ + 1\ (note that (h n j, h n j) = \h n j\ 2 = 1). □ 
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Lemma 2.4. The collection of asymptotics 



N 



g-* = 2vrV(l + ^) for all j = 1,2,.., N and ^ P nJ = I N + i\ 

3=1 

is equivalent to 



B n = 

Proof. As in Lemma 12.31 we set H n 
Note that B n = H n H* while 



2ttV(/ + ^). 



1 

2 h 



) 9n,N n n,N 



H*H n = [g n ]g n t ■ (K ik , h nJ ) 



N 



j,k=l 



Thus, as above, asymptotics = 2ft 2 n 2 (I N + ^) and = 2n 2 n 2 {I N + ^f) are 

equivalent. The diagonal entries of H*H n are g~j, so gr~J = 2ft 2 n 2 (l + £ 2 ). Asymptotics 



of the non-diagonal entries give (h n ^,h n j) = 2K 2 n 2 g 1 J l 2 :g 1 J l \ ■ i\ = i\, j ^ k, which is 



equivalent to J2j=i Bn,j — In + ^\ due to Lemma 12.31 
Note that, for sufficiently large n, 



□ 



A? 



A— A„ - 



2m 



M(X)d\. 



|A-7r 2 n 2 |=3|| 



This formula allows us to determine sharp asymptotics of B n (V). Moreover, it defines 
the analytic continuation of B n (V) for non-self adjoint potentials. 

Proposition 2.5. The following asymptotics hold true 

BJV) = 2ft 2 n 2 I N - — {{l^tjV}^ + O 

uniformly on bounded subsets of potentials V G £ 2 ([0, 1]; C NxN 
Proof. It's well known that 

x(0,z 2 ,V) = ^(l,z 2 ,V i ) 



+ 



^^/iv + i/ sinz(l-t) • V\t)sinztdt 



I dxsinz(l-x) ■ V\x) [ smz(x-t) ■ V\t) sinztdt + o(- 
2 Jo Jo V 



I Imz\ 



uniformly on bounded subsets of V. Substituting z = ft n + /x, = 3||V|| = 0(1 



one obtains 

-2„2 



x(o,^+»,v)= { -JT£i N + 1 



+ 



+ 



2ft 2 n 2 
-i 



7T 2 n 2 



7T 

+ 



sin 7m(l —t) sin imt ■ V\t)dt 

/ ((1— t) cos7rn(l — t) sinrmt + tsmftn(l — t) cos ftnt) ■ V\t)dt 
Jo 

— [ dx [ sinnn(l — x) sin nn(x—t) sin rrnt ■ VHx)VHt)dt 
n A Jo Jo 

1) 



2ft n 



n' 



2ft 2 n 2 



l±K n + L n + 
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where the matrices 



Znn lixn 



L n = -y»(°) + vi {cn) + o - = -y (0) + v {cn) + o I 

do not depend on /i. Hence, if /i = 311^11 and n is sufficiently large, then 

[ X (0, 7T V + /i, V)]- 1 = [/i^ n + L,,]" 1 + O ( 1 



27r% 2 



Also, note that 



1 /-i / |Im«[\ 

x '(0,2 2 ,\/) = -^(l,* 2 , V*) = -coazI N --J cos z(l-t)-V\t) sin ztdt+O I J 



Therefore, 



X '(0,7rV+/i,V) = (-l) 



n-l 



2 7m \ n J 



and 



2vr 2 n 2 



'AT 



2ixn 



y(sn) 



K?\pI N + I*K?]- x + o(^j. 



Since L n K n x doesn't depend on /x and 3||V|| = > \L n K n x \ for sufficiently large n, 
we have 



1 



n - 1 



— (f) [fiI N + L n K n x ] dfi = I N , 



M=mv\\ 



and so 
1 



2ir 2 n 2 n 



in 



2-nn 



y(sn) 



^ 1 + o[^)=iN- = [{i-t)v\™ + o(±). □ 



9 I -i> LV / J ' \ 9 



2.3. Proof of the direct part in Theorem 11.11 

Proof. In fact, all needed asymptotics have been obtained in Sect. 12.11 12.21 First, 
asymptotics of the eigenvalues and the individual projectors have been derived in 
Proposition 12.11 Second, asymptotics of the norming constants and the averaged pro- 
jectors follows from Proposition 12.51 and Lemma 12.41 In order to prove (C) suppose 
that £ : C — > is some entire vector- valued function such that P a ^(X a ) = for all 
a^l, £(A)=0(el ImVX l) as |A| -> oo and £ G £ 2 (R+). Due to Lemma 2.2 |CK06b| . 

[ X (0A,V)]- 1 = [<p*(lXV)]- 1 = (Z- 1 + 0(\-\ a ))((\-\ a y 1 P a + P^) as A — > A Q 

for some Z a such that det Z a ^ 0. Hence, the (vector- valued) function 

a;(A) = [ X (0,A,\/)]- 1 £(A) 

is entire. It follows from (12.11) that 



a; (A) = 0(|A| 1/2 ) as |A| = ir 2 (n + ±) 2 



oo. 



Thus, the Liouville Theorem gives u)(\) = u)(0) = uj q G and £(A) = x(0, A, V)cu . 
If co»o 7^ 0, then this contradicts to £ G £ 2 (M + ) in view of asymptotics (12.11) . □ 
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2.4. Explicit formula for the Weyl-Titchmarsh function. In this Sect, we prove 
that the Weyl-Titchmarsh function M(A, V) can be written as the regularized sum over 
all its poles. In other words, we give the explicit formula for M(A, V) involving only 
the spectral data X a (V) and B a {V) = — ies \ = \ a M(\, V). The proof is quite standard. 

Proposition 2.6. Let V = V* E C 2 

N 

E 

E a 



,l];C NxN ) satisfy ( TOj) . Then 



M(A) 



A-^coty^-^ 



27r 2 n 2 P° 



A„ — A 



N r 



+00 

EE 

n=n° j'=l 



5 



27r 2 n 2 P° 

TT 2 n 2 + Vj—\ 



(2.4) 



L a=l " n=l j=l ' 3 

T/ie series converge uniformly on compact subsets of C that do not contain poles. 
Proof. Note that 



Ay (A) 



B 



a -j 



A n ,j — A 
5 



27r 2 n 2 P° 
TT 2 n 2 +Vj — A 
2vr 2 n 2 P° 



2vr 2 n 2 P ? °) 



2^2 



-7T n 



n 2 n 2 



A 



(vr 2 n 2 - A) (vr 2 n 2 + vj - A) (A nj - - A) (7r 2 n 2 + v° - A) 



Due to Proposition 12.51 for the first terms one has 



N 



B ,_ 2vr 2„2po d -2n 2 n 2 I 

Z — ' 7v n — A 

i=i 



iV 



ir 2 n 2 



where (i n )^ n <> G ^ 2 . In particular, the series 
side singularities. Moreover, 

+oo +oo 



+oo 



-A 7T 2 n 2 — A' 

Pn^(A) uniformly converges out- 



E 4"(A) < 1 E 



n - 



as |A| = 7r 2 (m+i) i 



oo. 



Since P n j = 2n 2 n 2 {P° + £ 2 ) and A r 



7r 2 n 2 



+ v°: + £ 2 , the similar results hold true for 



'11,3 ~" " v- 1 j ' " y «j ■■ ■- ■ 

the sums of second and third terms of D n j(X). 

Thus, the right-hand side of ( 12. 4ft converges outside singularities and tends to zero 
as |A| = 7r 2 n 2 (m+|) 2 — > oo. It follows from the standard asymptotics of fundamental 
solutions that the left-hand side of (}2.4j) also tends to zero as |A| = -n 2 n 2 {m+^) 2 — > oo. 
Since the residues of both sides at singularities coincide, ( 12.41) holds true for all A. □ 

3. Inverse problem 

3.1. Proof of the surjection part in Theorem 11.11 General strategy. 
Step 1. Let some data (A* , pt, gjja^i satisfy conditions (A)-(C) in Theorem II. II and 
Bit = Pa(.9a)~ 1 Pa ( we use different superscript o for eigenvalues in order to make the 
further presentation more clear). Consider eigenvalues A* (possibly multiple for several 
first a). One can split them into N simple series {A* ,-}^Li, j = 1, 2, .., N such that 

{K,,l}n=l U {KiAn^l U U {A* iN } n =i = {Kx}a^l 



(counting with multiplicities) and A*j 



ln=l 

vr 2 n 2 + v ° + I 2 for all j 



1,2, 



N. 
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Using the well known scalar inverse theory (see fll.4p ) we construct some scalar 
potentials u|- G £ 2 ([0, 1]) such that 

f v^(t)dt = v° and a(v%) = {K,M=i ■ 
Jo 

Note that the corresponding isospectral sets are infinite dimensional manifolds, so there 
are infinitely many choices for each Vjj. For technical reasons, we choose i>*- such that 

9n l i. v< jj) = ~ res m (^ v jj) = 2vr 2 n 2 for all sufficiently large n, 

where m(X,Vjj) is the Weyl-Titchmarsh function of the scalar potential and 

X '(0,A a ,4) ^0, i.e., m(A Q ,^.)^0, for all a ^ 1. 

(one can always choose such i>J. in two steps: taking the scalar m-function with all 
residues equal to —2ir 2 n 2 and changing the first residue slightly in order to guarantee 
m(X a ,Vjj) 7^ for all a ^ 1). Let 

Thus, cr(V ) = {\a} a ^i counting with multiplicities. Denote 

Bi=Pi{gir\pir = B a {v«). 

Since V° is a diagonal potential, each subspace is spanned by some (one, if a is 
large enough) standard coordinate vectors e° and all are coordinate projectors. 

Step 2. Let 

A a (v) = M-\xi) = [ x (xT 1 mK,v) 

and 

Ai 1 = plA a (piy : El - C A' 2 = plAM)* : - C 

^ = <£A*(p*J* = £ - ^ 2 = <f a A a {q« a y : - &, 

where : — > 5*, : — > (E^) 1 " are the coordinate projectors. Note that 
A 1 a 1 (V ) = O, A 12 (V°) = 0, A 2 a 1 {V <> ) = and det A 22 {y ) ^ for all a ^ 1 

due to plx(0, O = AJ, n(P°JT = and det X '(0, A», V*) ^ 0. 

In order to describe some neighborhood of the isospectral set Iso(V°) near V , we 
introduce k a x k a matrices (more accurate, operators in the coordinate subspaces £*) 

A a (V) = K - A 12 (A 22 )~ l AH] (V), (3.2) 
Then (see Proposition 13.21 and Lemma [3. 3 p 

(i) all A a (V) are well-defined in some complex neighborhood B(V°,r°) of V°; 

(ii) for V = V* E B(V\r°) one has A a {V) = [A a (V)]* and the following holds: 

A a {V) = iff A* an eigenvalue of V of multiplicity k a . 
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Furthermore, for potentials V sufficiently close to V°, we set 

B a (V) = -^-(f M(X,V)dX, where cf = |min Q ^i(A* +1 -A*) > 0. (3.3) 

Zlxl J\\-\%\=d° 

If k° = 1, then M(A) has exactly one simple pole inside this contour, so B a (V) =B a (V). 
If > 1, we do not know precisely how the multiple eigenvalue A* is split, so B a (V) 
denotes the sum of all corresponding residues. Then (see Proposition 13.21 Lemma 13.31) 

(i) all B a (V) are well-defined in some complex neighborhood ^(V^r ) of V°] 

(ii) for V = V* e B(V°, r°) one has B a = B* a , rank B a = and the following holds: 

A a (V)=0 B a (V) = B a (V). 

In other words, B a (V) is the analytic continuation of B a (V) from the isospectral set 
Iso(V°) into some complex neighborhood of V° (emphasize that, due to the possible 
splitting of the eigenvalue A* in case fc* > 1, the original function B a (V) is discontin- 
uous even for self-adjoint potentials close to V°). 

Step 3. We introduce the mapping 

$:V^(A a (V);B a (V)) a> i 

which is defined in some complex neighborhood BiV ,^) of V° (see Sect. 13.21) . We 
prove that $ maps i3(V r<> ,r <> ) into some "proper" £ 2 -type space. In order to have the 
"nice" description of the image space, we consider some modification $, see details in 
Sect. 13.31 13.41 The modified mapping $ is analytic in BiV® , r°), so its restriction onto 
self-adjoint potentials close to V° is real-analytic. Note that, if V — V*, then both 
fc*xfc* matrix A a and NxN matrix B a , ranki? a = are self-adjoint. So, the total 
number of (real) parameters in (A a (V), B a (V)) is {k° a ) 2 + A£(2iV-A£) = 2Nk%. 

Step 4 We check that the Frechet derivative dy>Q of the modified mapping $ at the 
point V° is invertible (see details in Sect. 13.51 13.6P . Therefore, due to the Implicit 
Function Theorem, for each sequence sufficiently close to (_E>^) Q ^i there exists 

some potential V (close to V°) such that A a (V) = Ajy*) = and B a (V) = B* a 
for all a ^ 1. If a' is large enough, then the sequence 

B' a := B° a , if a ^ a*, and B' a := if a > a*, 

is close to (B^) a ^i. Thus, we obtain some potential V* such that 

A a (V')=0 for all a ^ 1, i.e., <r(V) = {\%} a>1 

(counting with multiplicities) and 

B Q (V) = B Q (V) = Bl, for a > a\ 

Finally, using the isospectral transforms constructed in |CK06bj . we change the finite 
number of residues B a , a = 1, 2, .., a' (see details in Sect. 13.71) . and obtain the potential 
having the given spectral data (A„,S^) a ^i or, equivalently, (A*, P^, g^) a ^i- D 



14 



DMITRY CHELKAK AND EVGENY KOROTYAEV 



3.2. Rough asymptotics of A a (V) and B a (V). This section contains some pre- 
liminary calculations. Loosely speaking, we consider the diagonal potential V as the 
unperturbed case and derive some rough asymptotics of spectral data for V close to V"° . 
The main results are formulated in Proposition 13.21 and Lemma 13.31 

Let <^°, ft®, x°j V° De the standard diagonal matrix-valued solutions (recall that 
V° is diagonal) of the equation —i])"(x) + V r<> (x)^(x) = \ip(x) satisfying the following 
boundary conditions: 

#>(0) = (^)'(O) = I N , rf(i) = -(x°)'(l) = Is, 

(r)'(o) = ^(o) = o, (Tf )'(i) = x °(i) = o. 

We denote ^Jx) — i P <> { x , Kt)i $a( x ) = ®{ x i Kt) an d so on - Let 

r( x ,t) = <p«(x)#«(t)-r(x)<p«(t) = - x o( x ) v o( t ) + v < >(x)x <> (t) 

be the (diagonal) solution of the same equation such that J°(t,t) = 0, (J°)' x (t,t) = Is- 
Let V = V° + W be some complex potential close to V°. Then x(x, A, V) can 
be easily constructed by iterations with the kernel J°(x,t) (note that |J°(x,t;z 2 )| = 
^~i e \imz\-\x-t\^ s t ar ting with j£(x, A) . Thus, 

X (0, z\ V) = X °(0, z 2 ) + J <p*(t, z 2 )W(t) X «(t, z 2 )dt + O [ L -j^ 3 J , (3-4) 

X '(0, z 2 , V) = ( X *)'(0, z 2 ) - J r(t, z 2 )W(t) X \t, z 2 )dt + O ( MLAJL J (3.5) 

uniformly on bounded subsets of W. In particular (see I2.2p . if \i = 0(1), then 
x(0,Kj+l*,V) = ^(diag{ yU -^+ Wl ,..,/i-^+<}+o(l)+0(||iy||)), 

X '(0, AV+/i, V) = {-l) n ~\I N + Oin- 1 )) , as n - oo, 

uniformly on bounded subsets of W. Recall that A a (V) = [x(x') -1 ](0, A* , V) and its 
block A 22 = qlA a {q%)* are given by (O) and d° = \ min a> i(A* +1 - A*) > 0. 

Lemma 3.1. There exists r° > such that for all (possibly non-self adjoint) potentials 
V e B{V°y) = {V e £ 2 ({Q,l};C NxN ) \\V-V°\\ < r°} 

the following is fulfilled for all a ^ 1 : 

det x'(0, A^, V) ^ 0, det A 22 (V) ^ and det x(0, A° +//, V) ^ 0, i/ = tf . 
Moreover, for all j = 1,2, ..,N and \fi\ = d° , 

[AZiV)]' 1 = 0{n 2 ) and [ X (0, X^+^V)]' 1 = 0(n 2 ) (3.7) 
uniformly on BiV® , r°). 

Proo/. It follows from ([311]) that all matrices x'(0, A^-, V), A 22 -(y), x(^,K,j+^ V ) 
are non-degenerate and ( 13 .TP holds, if n ^ is sufficiently large and r° is sufficiently 
small. So, one needs to consider only some finite number of first indices a = 1, 2, .., a*. 

Note that det X '(0, \° a , V ) ^ 0, det A 22 (V°) ^ 0, det X (0, A*+/i, 7*) ^ for all a 
and all these matrices (as functions of V) are continuous at V°. Therefore, if ||W|| ^ r° 
and r° > is small enough, then all x'(0, A*, V), A 22 (V), x(0, A* +//, V), a = 1, 2, .., a*, 
are non-degenerate too. □ 
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Proposition 3.2. (i) There exists r° > such that all A a (V), B a (V), a ^ 1, are 
well-defined by ( TO)) . ( TO) and analytic in i3(^°,r <> ). 

(wj For = 1,2, ..,/V iae asymptotics 

A n>j (V) = o(^iy B n>j {V)-B^ = 0(n 2 e n {W)), s n (W) = \W^\ + 
hold true uniformly for potentials 

V e B ^,^) = iv = V°+W e B{V°y) : / W(t)dt = 



Proof, (i) Due to Lemma I3.1[ all A a (V), B a (V) are well-defined in some complex 
neighborhood BiV^^r^) of V . These functions are analytic in this neighborhood since 
x(0, A, V) and \ V) are analytic for each A as functions of V. 
(ii) Let A = 7i 2 n 2 +[i and \n\ = 0(1), thus 

<f°(t, A) = (rrn)- 1 sin Tint + 0(n~ 2 ) and (-1)"~ V(*, A) = (7m) -1 sin 7mt + 0(n~ 2 ). 
Using (E3D, (ESD and £ W{t)dt = 0, we get 

X (0, A, V) = X °(0, A) + O (^P) , x'(0, A, V) = ( X °)'(0, A) + O f H 

(note that n -1 ^!! ^ e„(W / ) by definition). Due to (13. 6p . it gives 



AnAV) = IxixTW, a: ,, v) = Ani(v*) + o 



n 2 



Since A^V)^, ^ 2 .(V r<> )=0, A^(V r °)=0 and (^.(V))" 1 = 0(n 2 ), we have 



Due to the similar arguments, if A = A£ ,•+//, |u| = d°, then 

[x'x^Ko, a, v) = [(*W) _1 ](o, A) + O (n 2 e„(iy)) . 

Integrating over the contour |u| = gP, we obtain B n j{V) = B° n - + O {n 2 e n (W)). □ 

Lemma 3.3. For some r° > and all V = V* E i3(V°, r ) iae following hold: 

(i) Aa(V) = [A a (V)}*, B a (V) = [B a (V)}* and mnkB a (V) = k%; 

(ii) A a (V) = if and only if A* is an eigenvalue of V of multiplicity k^; 
(in) ifA a (V) = 0, then B a (V) = B a (V). 

Proof, (i) If V = V*, then M(A) = |M(A)]*, A G C. In particular, B«(V) = [B a (V)]*, 
Ajy) = [A a {V)]* and A a (V) = [A a (V)\*. Due to Lemma EH detx{0,X,V) has no 
zeros on the circle |A — A a | = d* for all V G BiV^^r^). Since the spectrum depends on 
the potentials continuously, for each self-adjoint potential V = V* G B(y°, r°) there 
are exactly fc* eigenvalues in the interval (A* — d°, A^+c? ) counting with multiplicities. 

If a > a , then fc* = 1 and rank5 Q ,(l / ) = iankB a (V) = 1. If a ^ a°, then 
rankB a (V) ^ Note that ranki? Q ,(l /<> ) = and i? a is a continuous function of V. 
Thus, if r° is small enough, then rank B a (V) ^ fc* for all a ^ a and V G ^(V , r°). 
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(ii) Recall that A* is an eigenvalue of V of multiplicity k? a iff dimKer x(0, A*, V) = k? a . 
Since det x'(0, V) ( see Lemma EH]) , this is equivalent to say that 

dimKer[x(x')~ 1 ](0, A*, V) = A£, i.e., mnkA a (V) = N — k a . 

Due to Lemma EH detAf(V) ^ for all V G <B(V*,r*)- Then, the last statement is 
equivalent to A a (V) = [A^ 1 - A 12 {A 22 )~ l Af}{V) = 0. 

(iii) If A a (V) = 0, then A* is an eigenvalue of multiplicity and there are no other 
eigenvalues in the disc |A — A*| < (P. Thus, 

B a (V) = -res A=A oM(A,y) = B a (V). □ 

3.3. Analyticity. Expanded mapping Proposition 13.21 (i) guarantees that all 
matrices A a (V), B a (V), a ^ 1, are well-defined in some neighborhood BCV®, r°) of V™. 
Let a° ^ and n° ^ 1 be such that 

k\ + kl + .. + kl« = N{rf-l) and A£ = 1 for aU a > a +l, 

so the double-indexing (n,j), j = 1,2, ..,N, is well-defined starting with n°. Also, let 
n° be sufficiently large such that <7nj(V°) = 2n 2 n 2 for all n ^ n° (see Step 1 Sect. 13. ip . 

Recall that £ n (V) = Ef=i^n,i(^) f o r « > 
Definition 3.4. Introduce the (formal) mapping 



= [A a -B a 



-oo 
,_„o I ) 



; 7m 



2 7 r 2 r i 2 

' j— j- 

(!) ,Tr( 2 ) 



Note that W a / and W n > map i3(V , r°) into some finite-dimensional spaces. Namely, 

^^(r.rVC^eC^ and *?:B(r,rVC N ©P x Y®C ra . 

Since has the finite number of components, it also acts into finite-dimensional 
Hilbert (Euclidian) space H {1) = 0^ [C fc « xfc « © c 7VxiV ] . It has been shown in Sect.ES 
that the components of \1/^ 2 ^ have "nice" asymptotics for potentials 

V e B°(FV°) = jv = G B(F°,r°) : jf = o|. 

Let N n » = {n6N:n^ n 9 } and C™ xm = {A = A*<E C mxm } be the real component of 
the complex Hilbert space C mxm , i.e., the real space of all self-adjoint mxm matrices. 



Lemma 3.5. (i) ^ maps B (F°,r°) irate ft( 2 ) = £ 2 C (^N n <> ; C^© [C JVxJV ] iV © C 
Moreover, the image [B°(V <> , r°)] bounded in TC^ 2 \ 

zs an analytic mapping between complex Hilbert 
spaces. Moreover, the Frechet derivative dyo'ty of ^> at V° is given by the Frechet 

derivatives of its components: {d v «^)W = (((d v ^ ( a ] )W)^ =1 ; ((4<>*l 2) )^+™o) . 
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(Hi) ¥ : B^V^r*) = B°{V\ r») n £ 2 ([0, 1]; C^ x7V ) ^ H R = W« x H$ is a real- 
analytic mapping between real Hilbert spaces and the Frechet derivative dy°^ is given 
by the Frechet derivatives of its components, where 

= ©1 [ c « Xk " © c m X 1 > = 4 ( ; ^ © [Cf xjv ] "© C 

Proof, (i) Due to Proposition 13.21 for all j = 1, 2, .., N 

A n>j (V) = 0(n- 2 e n (W)) and B nJ (V) - B* d = 0(n 2 e n (W)) 
uniformly on BiV® , r°), where 

e n (W) = \W^\ + JE1 so V] knW| 2 = 0(\\W\\ 2 ). 

n=n° 

Since B° nj = (cfr^Pf = 2ix 2 n 2 P^ n ^ n\ we obtain 

Wi nj (y)) +M G£ 2 and f^#-^ y°° Gf, j = l,2,..,N. 



n=n" 



uniformly on i3°(V°,r ). Also, due to Proposition 12.51 

B (V) l\ +oc 
7rn -^f -J* G £ 2 uniformly on S (F o ,r°). 

L 2vr% 2 J ; n=n0 

(ii) Due to Proposition ^. 21 all coordinates , a = 1, 2, .., a , are analytic in £>(V° ', r°). 
Hence, yfr^ is analytic too. Similarly, all coordinates n ^ n°, are analytic in 

BiV®, r°). It follows from (i), that ^Z 1 - 2 -' is also locally bounded in B°(V°, r°). Therefore 
(e.g., see [PT87j (Appendix A, Theorem 3) or [Di99| (Chapter 3, Proposition 3.7)), 
\l/( 2 ) is analytic as the mapping between Hilbert spaces and its Frechet derivative (or, 
equivalently, gradient) is given by the Frechet derivatives (gradients) of its components. 

(hi) By Lemma \I/ maps £>^(K , r ) into 7Yr. \I/ is real-analytic due to (ii). □ 

3.4. Analyticity. Modified mapping <fr. The expanded mapping \l/ introduced in 
Definition 13.41 is real-analytic but overdetermined. In other words, its coordinates, 
obviously, are not independent from each other. In particular, there are no chances 
that the Frechet derivative dyoty is invertible. On the other hand, the coordinates 
A a (V), B a (V), a^l, of the original mapping $ are independent, but we have no 
"nice" description of the image space. The next goal is to construct some modified 
mapping $ = ($( 1 ),$( 2 )) (see Definitions I3T61 f378l I3T91) such that 

(i) it keeps the full information about A a (V), B a (V), a^l; 

(ii) it is real-analytic as the mapping between Hilbert spaces; 

(iii) its coordinates are "independent" from each other (more precisely, 
in Sect. 13. 5\ I3.6l we will show that dy©$ is an invertible linear operator). 

We start with a slight modification of the first coordinates B a (V), a = 1, 2, .., of. 
Recall that, if V e B^(y*,r°), then B a (V) = [B a (V)]*, vankB a (V) = k» and 

B" a = (plTBlpl, pIBXpD* = {glY 1 = [{gl)' 1 }* > 
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(moreover, g° a is diagonal, since V is diagonal). Therefore, if r° >0 is sufficiently small, 
then for each a = 1, 2, .., a we have the (unique) factorization 

B a = MT + (<£)*£«] C a \pl + Kql] , ° a = ° a u Z*\ ' S ° ~* f"' , (3.8) 

where B l a l = p%B a {p%y, Bf = q°B a (pl)* etc. Note that C a > 0, since rank5 a = k%. 
Definition 3.6. We introduce the first component of the mapping $ by 

: Bl (V, r°) - wj? = 0tJc^ x ^©C^ xfc2 ©C (7V ^ )x ^, 

(3 9) 

$(D(V) = ($«(V))^ =1 , = (A a {V) ; C a (V) ; E a (V)) . 

Remark 3.7. Z>ue to Lemma l375\ (ii), $^ is well-defined and real-analytic in ^(V°, r°), 
z/ r° > z's small enough. Note that $W can 6e reconstructed from $W and tee tote/ 
number of real parameters containing in $W is 2N{k\ + k% + .. + k^<>) = 2N 2 (n° — l). 

We pass to the design of the second component $( 2 ). The main purpose of (rather 
technical) Definition [331 is to combine heterogeneous objects from (11.61) into one object 
having "nice" asymptotics as n — > oo (see Proposition 13. 101) . 

Due to Proposition 13.21 if r° > is sufficiently small, then 

\KA V )\ = 0(n- 2 e n (W)) and \B nJ (V) - 2vr 2 n 2 P J °| = 0(n 2 e n (W)). (3.10) 

In particular, if V G ^(V^ , r°), then factorization (13.81) is well-defined for all n ^ n°. 
Recall that fc* ■ = 1, so A n j(l/) and C n j(V) > are real numbers. 

Definition 3.8. Let V G B^V®, r°) and r° > &e sufficiently small. Introduce two 
numbers a n j(V),c nt j(V) G K and one vector e n j(V) G sttcn teat (e n j,e°) = 1 as 

a nJ (F) = 2vrVl nj (\/), c^fV) = [(27r 2 n 2 )- l C n!j (V)]\ e n>j (V) = e° + K J (^)e J °. 

Furthermore, define NxN matrix Y n = Y n (V) G C A?xAr by 

Y n = ( exp[te n> i] ■ Cn,i ■ e n ,i ; exp[ia„ i2 ] • c„ i2 • e n>2 ; ••• ; exp[te„ iA r] • c UtN ■ e n , N ) 
and let 

Y n (V) = U n {V)S n (V), U* = U~\ S* n = S n > 0, 
be its polar decomposition. 

Note that all A n j, B n j, j = 1,2, ..,N, can be easily reconstructed from U n , S n . 
Factorization (13.81) reads now as 



Icy 2 2\-ln _ 2 * 

{Z7i n ) n n j — c n j ■ e n ,j e n ,ji 

so (13. lOf) gives 

K,j(V)\ , \cnj(Y) - 1| , |e nJ (y) - e°| = 0( £ „(W)) 
uniformly for n ^ n°. Hence, 

|y„(v) - i N \ , |^ n (y) - iiv| , \s n (v) - i N \ = o(e n (w)) (3.11) 

uniformly for n ^ n° and det i^(V) 7^ for all G i3^(V^, r°), if r° is small enough. 
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Definition 3.9. Formally introduce the second component of the mapping $ by 

$(2) . V ^$(2)(^/) = ($f (V))+- , 

$( 2 ) = (-HogU n ; 27rn-(S n -I N )) : B°(V, r<>) - C^ xiV © C^ xAr , 
w/iere logC/ n = (U n -I N ) - \{U n -I N ) 2 + §(£4-//v) 3 - ... 
Recall that A nJ (V°) = and B n , 3 {V°) = 2n 2 n 2 P° for all n ^ n°. Thus, 

^(K ) = U n (V ) = S n (V ) = I N and ^{V ) = (0 ; 0) for all n ^ n°. 
Proposition 3.10. There exists r° > stzc/i t/iat £/ie mapping 

^ : r«) - 4(N B o; N x 

zs well-defined and real-analytic in £^(7°, r°). Moreover, the Frechet derivative dv>& 2 ' 
of $( 2 ) at 7* is given fry £/ie Frechet derivatives of its components. 

Proof. Due to fl3TTTD and |e„(17)| 2 = 0(||jy|| 2 ), for sufficiently small r° > the 

mapping 

y-.V^ (Y n (V)-I N )+~>, B R (V»,n - 4(N n .;C"*"), 

is well-defined. Recall that Y n is some simple function of A n j and B n j, j = 1,2, ..,N 
(see Definition 13.81) . Using real-analyticity of the first two components of the expanded 
mapping (see Definition 13.41 and Lemma [3. 5p . we conclude that y is real-analytic 
as a composition of real-analytic mappings. Since S n = (Y*Y n Y^ 2 and U n = Y^S 1 " 1 , 
both mappings 

and 

u : u _> (-ii og c/ n (y))+-o, B£(VV) - 4(N„.;C?*"), 

are real-analytic too as compositions of y with some simple coordinate-wise transforms. 

In order to complete the proof it is sufficient to show that S actually acts into 
"better" space £ 2 . Note that 

1 A= B n 



Y n Y n - ^2 c n j ■ e n je n j - B 
j=i 3=1 

Due to Lemma 13.5} the mapping 

Z : V - 2vm • (F„F; - 1*)+^, BjtfV, O - £(N n o, C^ xAr ) 

(which is the third component of \]/( 2 )) is real-analytic. Using S n = [U^ 1 (Y^Y*)!!^ 1 ' 2 , 
we obtain that the mapping 

S:V^2nn- (S n (V)-I N )+Z», B&(W*) - £(N n .; C^ xjv ), 

is real-analytic as a result of some coordinate- wise transforms with Z and U. Note that 
$( 2 ) = j S). Since the Frechet derivative dyo^ is given by the Frechet derivatives 
of its components, the same holds true for all mappings y, S, U, Z and S. □ 
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Remark 3.11. The mapping $ = ■ <|>( 2 )) i s real- analytic too, since both $( 2 ) 
are real- analytic, and its Frechet derivative is given by the Frechet derivatives of 

n n°, contains 2N 2 real parameters, i.e., exactly 
"the same amount of information" as, say, the n-th Fourier coefficient V^ n ' . 

3.5. Explicit form of the Frechet derivative d v «<3?. We denote by 

V° : W{x) i — ► W(a;) - JF (0) 

the orthogonal projector in £ 2 {[0, l];C% xN ) onto {IF G £ 2 ([0, 1]; C^ xiV ) : = 0}. 

Recall that the mapping $ was introduced in Definitions 13.61 and 13.91 Due to 
Remark ETH (d V o$)W for W G P°£ 2 ([0, 1]; C^ x7V ) is given by 



a 

n > n°. 



1,2,.., a* 



(dyA a )W, (rfyoC Q )W, {d v »E a )W for 

and (d V o[/ n )W, (dv»S n )TF for 

We need some preliminary calculations. Let 

Xt = xi-iKiV ), <P% = <p(;X%,V*) and so on. 

Since V° is a diagonal potential, all these matrix-valued functions are diagonal. For 
short, we will use (a bit careless) notations like 

:= xiimxino)]- 1 = [(xirm^xm. 

Recall that p^ : — > £* and : — > (i^) -1 are some coordinate projectors. 
Note that Ker[£(0)(£)*] = {0}, Ker[( X :)'(0)(p»)*] = {0} and Ker|ft(0)(p*)*] = {0}. 
Thus, expressions 

[^(0)]-W, [(Xiy(0)}-\PIT and [^(O)]- 1 ^)* 
(and their conjugates) are well-defined. 

Proposition 3.12. For all a ^ 1 and IF G P°£ 2 ([0, 1]; C^ xN ) the following hold: 



and 



(d v «c a )w = p; 



(d v .A a )W = p, 
{d v ,E a )W = -q, 



X ° a ^ Wit) dt 
Wit) Xa ' dt 



1 xl(t) 



(p°J* 



W{t)^\ + W{t)^ 



where 



&(0) 



&(0) &(o) 

X°J0) 



(Pa) 



2x°a(0) 



*«(*)■ 



(3.12) 
(3.13) 

(3.14) 
(3.15) 



Proof. It follows from (13.4)) and (13.51) that 



(aW(0X))W= / ^(t)W(t);£(t)dfc, (W(0,A*))W = 



i?*(t)W(t)x«(t)dt 



and 



(d v «x(o,K))w 



mt)w(t) X i(t) + vi(t)w(t)r a (t)) dt. 
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Recall that A a = A^ 1 - A™ (A™)- 1 A 2 a l , where 

A Q (V) = ix(xr 1 ](0,K,V), A^=plA a (pir, A^=plA Q (qir and so on. 
Due to ^ 2 (U°) = 0, A 21 ^) = and p^xU ) = °> one obtains 

(d v «A a )W = (dvoA^)W = pl(d v » X (0, K)) W [(xl)\0)r\piy 



Pi 



(pi 



This gives (EH, since = p*x«(*)[(Xa)' '(a)]" 1 ■ Next > 

1 



(d v «B a )W 



2m 
1 

2ttz 



|A-A£[=d<> 
|A-A*|=ct 



(dv^xx-'mvwdx 

-(d V o X '(0, X))W + (dvx(0, A))W 



dA 



X°(0,A) ^" ^ ' " Jx»(0,A) 

Note that the diagonal matrix-valued function [x°(0, A)] -1 has the unique pole 
(at A^,) inside of the contour of integration and 

In nof In Xa(°) 1 no , ™ ^ 



K + Qt 



Ql + 0(X-Xl) asA^A c 



x o(o,x) a lxl(o)(x-x» a ) 2[xm} 2 \ a a x«(o) 

where Q° a = = - Recall that £ a = B^B^)- 1 and 5 21 (U* 

(dy.£ a )W = (d V oB 2 a l )W ■ [Bl'iV*)]- 1 = 

and 



0. Thus, 



(dvoB^W-pl-^-ipir 



(x«)'(o) 



(dv.B^W = q, 



x%(0) 



&(o) 



(Pa- 



using X«(0)^(t) + (^)'(OR(t) = Xl(t), one obtains (EZHJ). 

Furthermore, C a (V) = B£(V) = p^B a (V)(p^)* . In contrast to (dyoB^W, we 
do not have cancellations of the singularities by the projectors, so one should find the 
residue at the second order pole A*. Straightforward calculations give 

In 



(d V oC Q )W 



Pc 



res p° Q 



-(<W(o, x))w + (x 2?^ } (<W(o, a))w 



K(t)w(t) x i(t) + 



+ 



X°(0,A) 



&(o) 
ft(o) 



X°(0,A) 

v«(*W*)x«(*) 



(Po 



2[x a (0)] 2 



in 



J ft(0) 



XS(0) 
Using the identities 



— (y° y (q v (twft) ) 

WaJ W 2[^(0)] 2 



*(*) + Wrf(*) + W*(*)-* (t) 



X°a(0) ^ ' X°a(0) 

and p:(x:)'(0)<(f) = one obtains (GUI). 



Xl(0) 



□ 
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Introduce the functions 
X^(t) = [Xl{t)\n = X(t, A a , t&) and £*(f) = = £(t, A£, t&), 

where is given by (13.151) . 

Corollary 3.13. Let a ^ 1 and 1(a) = {s : A* G c(v^ s )} (by definition, the set 1(a) 
consists ofkl indices). Then, for all W G V°C 2 ([0, 1]; C^ xN ), 

[(d v ,A a )W] ]k = (W jk , ug fc >), [(d v ,C a )W] 3k = (W jk , j, k G 1(a), 

where for all A^ G cr(vjj) fl cr(v kk ) the functions u£ and Ua are given by 

u ( J k \t) = [(x^) / (o)(x° a fc ) , (o)]- 1 • x% j (t)x% k (t), 

(S.lb) 

u% k) (t) = [x^(o)x°/(o)]- 1 ■ {&(t)x°/(t) + xy(t)e a k (t)) ■ 

Furthermore, 

[(d v ,E a )W] jk = (W jk , uf>, j i 1(a), k G 1(a), 
where for all A* G cr(v kk ) \ cr(vjj) the function Ua is given by 

u { f } (t) = -ix^mx^rm- 1 . (3.17) 

Proof. Since x„> are diagonal matrices, this is exactly the result of Proposition 13.121 
rewritten in the coordinate form. □ 

Proposition 3.14. Let n ^ n° and j, k = 1, 2, .., N be such that j ^ k. Then for all 
W G V°C 2 ([0, 1]; C* xN ) the following identities hold: 

[(d v ,Y n )W] n = (47r 2 n 2 )- 1 (iy j , , 5<g>) + i • 2ttV<W0 , u^), (3.18) 

where the functions u^j and are given by ^3. Id) , and 

[(d v «Y n )W] jk =(W ]k ,uf k ) ), (3.19) 

where the functions u^ k are given by \3.1J) . Furthermore, 



(3.20) 



(d V oS n )W = \ ((d v .Y n )W + {(d V oY n )W}*) , 

(d v «U n )W = \ ((d v .Y n )W - [(d V oY n )W]*) . 

Proof. By definition of Y n , 

[(d V oY n )W} jk = ((d V o[exp(ia ntk ) ■ c Ujk ■ e n:k \)W , e°) . 

Recall that a n , k (V°) = 0, c^V*) = 1, e n , k (V) = e° k + E n , k (V)e° k and E n , k (V°) = 0. 
Thus, 

= (r\,r«r„ AW + 7: ■ (dr, A n- AW = ^ 

4vT 2 n 2 
and 

[(dyYnWj* = [(d v ^, fc )W1i • 
Due to Corollary ESI one obtains fl3~T8l and fl3~T9l . Recall that S n = (Y*Y n fl 2 , 
U n = Y n S~ l and Y n (V°) = U n {V) = S n (V°) = I N . This immediately gives (13301) . □ 



[(d v «Y n )W\ j:i = (d v .c nJ )W + i ■ (d v ,a nJ )W = ( dv * C ^ W + i . 2-K 2 n 2 (d v .A n . j )W 
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3.6. Invertibility of the Prechet derivative d v «3>. Due to Remark [3. Ill 
{dv&PW = {{d V oA a )W ; (d V oC a )W ; (d V oE a )W), a = 1,2,.., a\ 
(dv.^n^W = (-i(d v <>U n )W ; 27m(d V oS n )W), n = n*,n + l,. 



Recall that Wkj = Wjk for all l^k^j^N. It immediately follows from Corollary 13. 131 
and Proposition 13. 14l that the entries of the components of (dyo§)W are 

(1) for all j = 1,2, ..,N (diagonal entries of (a) A a , C a and (b) U n , S n ): 

(a) (Wjj , Ua ), (Wjj , Ua), where a < a are such that A* G o-(v^); 

(b) 2ix 2 n 2 ■ (Wjj , (27m)- 1 • (W n , Sgf> ; /or all n ^ n«; 

(2) for all l^k<j^N (non-diagonal entries of (a) A a , C a ; (b) E a ; (c) U n , S n ): 

(a) (Wjfc , Ua ), (Wjk,Ua ) and their complex- conjugates 
(Wjk, Ua 3) ), (W~~ k , Ua 3) ), where a < «°: A* G fl a(v£ fc ); 

(b) (VFj-fc , vf£ ), where a ^ a are snca i/iai A* G cr(^ fc ) \ a(Vjj); 

,(kj) 



(Wjk , u a ), where « ^ a* are such that A^ G c(^J,-) \ ^(^fcfc) 
.(?*) „.(*i)i\ _ /w/ r„,0' fc ) , „,( fc i) 



e,s 



(c) ^ ■ (W^ fc , KV-<j j ]), 7to ■ (W^-fe , [<V +<j ]) ana 1 taezr conjugat 
h-<Wk, [«£?-«£!?]>, ™ ■ <T^, [ujf +«g*J]>, /or a// „ £ n». 

Note that uL jfc) = and = v£ j) , if A Q G <r(u?.) n o(r^). 

Definition 3.15. For each 1 ^ k ^ j ^ N we introduce the collection of real scalar 
functions 



U (jj) = a ^ a o . A^G(t(4)} U {27rVng ) , (27m)- 1 u^, r^n 



U {u^\ a^a«:\%ea{vt k )\a(v« 3 )} U a<o? : A* G<r(^.) \ <r(t&) 

where the functions u£ and ni^' are given by A3. 16]) and {3.11 ). Note that each 
collection contains exactly 2(n°— 1) functions with "small" indices a a°. 

Remark 3.16. I>ue to tae arguments given above, in order to prove that [rfyo^] -1 is 
bounded, it is sufficient to prove that each V°U^ is a Riesz basis ofV°£ 2 (0, 1). 

Lemma 3.17. For each 1 ^ k ^ j ^ iV taere exists some collection of functions 
yO'fc) c pO£2( Q) ^ w ^ ic ^ - s ^orthogonal to U ijk) (and, therefore, to V°U ijk) ). 

Proof. Taking into account definitions ( 13. 16ft , ( 13. 17ft and ( 13. 15ft . it is sufficient to con- 
struct some collection yw fc ) C V°£ 2 (0, 1) which is biorthogonal to V°U^ k \ where 

= {j&jtf, for all A- G a(v»jj) U a{v% k )} 

U {r a J X° a k + X^xT, for all \» a G o(r^) n , 
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since U {jk) and U {jk) are related by some simple linear transformations (namely, multi- 
plications by fixed constants, (X)£ = X+ c x) (X; x) an d (^1,^2) (^1+^2, u\— U2)). 
Note that we consider both cases k — j and k < j simultaneously. Let 



o,j o,A*i/ 



r 0,1 <>■ 



for all A^ G cr(v?.0 U a(w 



-A; A:; 



U { W<pf + tfffi, for all \* G fl <t(<4) } 



By definition, V*C P°£ 2 (0, 1). Let \ a ^ \ p and {x, f} = XV* ~ xV- The standard 
trick (e.g., see |PT87j pp. 44-45 for the similar calculation in the scalar case ) shows 



(x?x5*)(^J*)' - {x^xT)'Wvf) (t)dt 



{X«a j ,^ 3 }(X°aV/) + {X°J^){XT^T} m 



o,k . _o,An 



(3.21) 



{x^^ j }{xii\vf}\\ _ k>;ki)-M(o) 



0,j 0,fc-| 



2(A a — A/3) 2(A a — \p) 

If both A Q , \ p G a{v%) U trKfc), then tf?{l)iptf{l) = X^^X^) = 0- Hence, 

x° Q 'V/, [^V/]') = 



Moreover, if A a G &(v"jj) H cr(t>£ fe ) (the case A/3 G H cr^/J is similar), then the 

right-hand side in (13.211) . as a function of A Q , has a double zero, so we can differentiate 
this identity (with respect to A Q ) and obtain 



(AT+x^r*, [$Vj*]') 



0. 



Also, if both A Q , A G cr(t£,) n a(vi 



then 



/ • 0.7 o.A 1 o,j • o,fe r -<> iJ 1 ^iJ • o.An/X 



Let X a = X p G ff(^)\<7(^) (or \ a = X p G \ a{v^)). Then 0&>>#} = 0, 

{X^,¥^'}^0and 



Let A a = A ( g G a«) R fffe). Then {x^,^} = {xr,<' fc } = and 



Using (I3.2ip for A/3 — > A Q , one gets 
Similarly, 



^'(1)^/(1) _ [y£WKl) 



lim , 

A/3^A Q 2(A a — A/3)' 



^0. 
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Finally, one needs to correct V^ fe ) slightly, replacing the functions [<fa'' ( Pa k + i P < a'P < a k \' 
for all X a G <t(v?j) D a{v° kk )) by 



with appropriate constants c Q , in order to guarantee 

After these corrections, V^ fc ) becomes biorthogonal to U^ k \ 

Proposition 3.18. V°U^ is a Riesz basis ofV°£ 2 {0, 1) for all l^k^j^N. 

Proof. Since V°U^ k ^ admits the biorthogonal system, it is sufficient to check that 
elements of V°U^ k ^ are asymptotically close (say, in £ 2 -sense) to some unperturbed 
Riesz basis (note that these functions are in one-to-one correspondence with eigenvalues 
of and v kk , and we have two functions in hi^ k ^ for common eigenvalues). Those 
u E that correspond to first eigenvalues A* -, A* k , n < n°, do not affect the 

asymptotical behavior, so it is sufficient to consider n ^ n°. 

We need some simple asymptotics. Let A = n 2 n 2 +fi, \x = 0(1), and v £ £ 2 (0, 1) 
be some (scalar) potential. Then 



□ 



. , . sin7rn(l— t) ^( 1 , 
x (t,\,v) = ^ '- + — , 



X '(0,A,t;) = (-ir 1 + O^ 
as n — > oo. In particular, 

t(t,\,v) = x(t,X,v) 

If k = j, one obtains 



x(0,A,u) 



(-D n 



;i-t)cos7rn(l-t) (\ 

2vr 2 n 2 U 3 ' ' 



2vr 2 n 2 



+ 



71" 



x(0,A,u) = O 



1 



X(0,A,,) x(tiA>v)= (l-t)cos7rn(l-t) +0 ^ ) 



2tt V • 



2*(0,A,u) 

2vr 2 n 2 [ X ^.(t)] 2 



[(x^)'(o)] 2 



2vr 2 n 2 



— cos 27rnt + O ( — 
n 



and 



(2~// ) J ;/ 



i . ~(n) 



-V° [(l-t) sin27mt] +0 



n 



It's easy to see that the collection 

TZ= |cos27rnt , V°[(l-t)sm2imt], n ^ l| (3.22) 

is a Riesz basis of V°C 2 (0, 1). Indeed, all functions (| — t) sin 27mt, n ^ 1, are linear 
combinations of cos27rmt, m ^ 1, since they are symmetric with respect to \. Hence, 

(/,cos27rnt)+ri \ _ ( I \ f (f,cos2Trnt)+™ 
(/,P o [(l-0sin27r^])+- ) ~ { A \I 



and the linear operator (/, cos 2imt) 
bounded in £ 2 , since the operator / i— > (h—t)f is bounded in £ 2 (0, 1). 



+oo 
n=l 



(/,sin27rnt)+ri 



(/, V [a-t)sm2nnt})^ f E £ 2 (0,1), is 
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Thus, 7Z is a Riesz basis of P°£ 2 (0, 1) and V°U^ j) is £ 2 -close to 1Z (note that 
in both r p°li^i> and 1Z there are exactly 2(n° — 1) functions with n < n°). Due to 
Lemma [3.1 7[ the elements of V U^^ are linearly independent. Therefore, V°U^^' is a 
Riesz basis of P°£ 2 (0, 1) by the Fredholm Alternative (see, e.g., |PT87j p. 163). 

Let k < j and n ^ n°. Due to [( x ^) , (£'4)- 1 ](0) = -(^ fc )- 1 = -27r 2 n 2 , one has 

Xnjfc (*) Xn%t) Xni (*) Xnl (* ) 



Note that 
and 



x&(o)(x#)'(o) 2n^xZmtM 

x;l(0) = (A; fe -A;,) • x°4(0) + 0(n" 4 ), 



since Xn,j(0) = anc ^ Xnj(O) = 0(n 4 ). Therefore, 



m 

n,k 



XnjWXnU*) 



+ 



+ - 



rr 



Thus, 



1 r (jfc) (fej')l 



and, since the first term of u^ k k \t) is antisymmetric with respect to j and k, 



cos 27mt „ / 1 
+ 

n 



A n,fc A n,j 



-P° [(l-t)sin27mt] + O 



r? 



As above, we see that V°U (jk) (up to some uniformly bounded multiplicative constants) 
is £ 2 -close to the Riesz basis TZ given by ( 13.221) . So, J>°U^ k ' is a Riesz basis due to the 
Fredholm Alternative and Lemma 13.171 □ 



Corollary 3.19. The Frechet derivative 

d v «$ = (d V o$ {1) ; d v «& 2) ) : P°£ 2 ([0, 1]; < 



lNxN^ 



n 



a) 



(2) 



(1) 



,(JV-A£)x*£ 



n 



(2) 



zs a linear isomorphism (in other words, dy>& is invertible). 
Proof. See Remark 13.161 and Proposition 13.181 



□ 



3.7. Completion of the proof. Changing of the finite number of first residues. 

Let {(A*, Pt, ^)}a>i be some data which satisfy conditions (A)-(C) in Theorem ll.ip 
and Pt = Pt(^t)-ipt. Reca n that P^ = P° + £ 2 and {g^y 1 = lTt 2 n 2 (\ + 
Similarly to Definition ^. 81 if n is sufficiently large, then we may introduce the (unique) 
factorization 



(27r 2 n 



Note that e. 



t . 

n,j 



e° + i 2 and d , = 1 + i 2 . Define 



"J 

C n,l ' e n,l i C n,2 ' e n,2 ) 



i C n,V ' e n,N 



iNxN 
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Since Y£ = In + £ 2 , the matrix Y^ is non-degenerate for all sufficiently large n, and we 
may introduce its (unique) polar decomposition 

Yl = ulsl plr = [u\]-\ [4]* = s!>o. 

Note that U ] n = I N + f and S\ = I N + P. By our assumptions, J2j=i p l,j = In + t\ 
and {gl^y 1 = 2n 2 n 2 (l +£f), so Lemma [23] gives 

uKsmuir = n<rtr = {^ 2 n 2 y i £f =1 = + <?■ 

Therefore, = J^ + f 2 . 

Recall that L£ = Ef n (V*) = J w and S° = S n (V°) = I N for all n ^ n* so $( 2 ) = 0. 
Since the Frechet derivative dyoQ is invertible, the mapping $ = $( 2 )) is a local 
bijection near l^ . Therefore, if a' is large enough, then there exists some potential 
V* G B^V ,^) such that 

= $i 2) (V) = $i 2) (^°) = for all n° < n < n% 

and $i 2) (^*) = (-*logt^ ; 2im ■ (Sl-I N )) for all n > n', 
where a'-a* = N(n* — rf) (i.e., a* + l corresponds to the double-index (n*, 1)). Since 
the original mapping $ can be reconstructed from $, one has 

l a (V) = A a (V°) = for all a < a*, 

2 nj .(V) = o and B nd (V) = B ] nj for all n > n*. 
Due to Lemma [3. 3 [ it gives 

<KO = KWi and S n j(V) = JB^- for all n ^n'. 

At last, we need to change the finite number of first residues (B a (V'))^L 1 to 
{Ba)a'=i- Recall that the isospectral transforms constructed in [CK06b] allow to mod- 
ify each particular residue B a in an almost arbitrary way. The only one restriction 
(concerning the change of projector P a to P a ) is 

.F a nRanP a = {0}, 

where JF Q , dimjF a = N — k a is some "forbidden" subspace that is uniquely determined 
by the spectrum and all other subspaces (£p)p^ a - It's not hard to conclude (see Propo- 
sition IA.4I) that this restriction is equivalent to the following: 

One can modify B a in an arbitrary way such that (C) holds true. 

In general situation one can change all B a (V) to B^ by a' steps. Nevertheless, it may 
happen that at some intermediate step the desired residue B^ violates (C). In order 
to overcome this difficulty note that one can always change B a to some B^ which is 
arbitrary close to B^ a in the natural topology. Then, in any case, after a' steps one can 
obtain some potential V* such that B a (V') = B^ for all a — 1, ..,a° (and, of course, 
B a {V') = B a (V) = Bl for all a > a 9 ). By Corollary \KM the set of all admitted 
by (C) sequences (B a )^ =1 is open in the natural topology. Therefore, if (B^)^' =1 and 
(Ba)a'=i are close enough, then all changes B^ h- > B^ are permitted. So, after another 
at most a* steps one obtains the potential V such that B a (V) — B^ for all a — 1, .., a' 
(and still B a (V) = B\ for all a > a'). The proof is finished. □ 
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A. Appendix. Property (C) 

Let A a > for all a ^ 1. Note that (C) doesn't depend on shifts of the spectrum, 
so we do not lose the generality. We begin with the following simple 

Remark A.l. If an entire function £ is bounded on the real positive half-line, then the 
condition £(A) = 0(e' Imv/ ^') is equivalent to say that i{z 2 ) is an entire function of 
exponential type no greater than 1 (see |Ko8 8j . p. 28). 

Recall that the Paley- Wiener space P^.y consists of all entire functions f(z) 
of exponential type no greater than 1 such that / G £ 2 (K). The Paley- Wiener theorem 
(see [Ko88] p. 30) claims 

fEPW^ iff f(z) = ^ / 4>(t)e- izt dt, where G £ 2 (-l, 1). (A.l) 

^ J-i 

Proof of Proposition [T73l If G £ 2 ([— 1, 1] ; C N ) is some vector-valued function 
such that 

J (j)(t)dt = and h* a j <p(t)e ±iVXZt dt = for all a ^ 1, (A.2) 



then 

— / <p(t)e~ lzt dt = zf(z) and P a f(±y/K t ) = for all a > 1, 

where zf(z) G PW hl>1] . Denote £(z 2 ) = l[f(z)+f(-z)} or ^ 2 ) = -/(-*)]. 
Then, P a £(A a ) = 0, a ^ 1, £(A) =0(el Im ^1) and £g£ 2 (M+). This contradicts to (C). 

Conversely, let f(A) =0( e l Im ^1) and £ G £ 2 (M+). Then /(z)=z£(z 2 ) G PW hl>1] , 
so it admits representation (lA.ljl with some G £ 2 (— 1, 1). It's easy to check that 
P a £{K) = and /(0) = imply flA~2l) . Hence, = 0. □ 

We have the immediate 

Corollary A.2. If one fixes the spectrum {X a }a^i and all projectors P a , a ^ a' + l, 
for some a' ^ ; then the set of all finite sequences (P a )^Li satisfying the condition 
(C) is open in the natural topology. 



Introduce the function 



X(0, A, V)Pl 

fr(A) = V (A.3) 

A — An 



where Pp : — > Si is the orthogonal projector onto the subspace Sp = Ker x(0, Xp, V). 

Proposition A.3. Let (3 > 1 and V = V* G £ 2 ([0, 1]; C* xN ). Then, 

(i) £p : C — ► C NxN is an entire matrix-valued function, £p(X) = 0(e^ lm ^) as |A| — ► oo ; 
G £ 2 (M+) and P a £p(X a ) = for all a^(3. 

(ii) If £ : C — > C N is an entire vector-valued function such that £(A) = 0(e' Imv/ ^') as 
|A| -> oo ; £ G £ 2 (M+) and P a £(A Q ) = /or a// a ^ /?, t/ien £ = £ p h for some h G C N . 
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Proof, (i) The function £p is entire due to x(0, Xp, V)Pp = 0. Furthermore, 

£g(A) = 0{\X\-h llm ^) as |A| -> oo and P Q ^(A Q ) = for all a ^ /3, 

since P aX (0,A a ) = P a [y? (l,A a )]* = [<p(l,\ a )P a ]* = 0. 
(ii) Lemma 2.2 [CK06bJ claims 

[ X (0 ) A,^)]- 1 = [^(l,A,V)]- 1 = (Z- 1 + O(A-A a ))((A-A Q )- 1 P Q + P Q ± ) as A - A Q 
for some Z a , a ^ j3, such that det Z a ^ and 

[ X (0, A, V)]- 1 = A, = ((A-A^P* + (P«) ± )(^ 1 + O(A-A^)) as A - A, 

for some Zp, det 7^ 0. Due to P a ^(A a ) = 0, a 7^ /3, the (vector- valued) function 

is analytic except Xp and u;(A) = (A — A/3)~ 1 P / §/i + 0(1) as A — > A/3 for some fa G C^. 
Since c<j(A) = OQAI 1 / 2 ) as |A| = 7r 2 (n + |) 2 — >■ 00, the Liouville theorem gives 

£(A) = X (0, A, V) [(A-A^Pj/i + cu ] = C/?(A)/i + x(0, A, V)u for some ^ G C". 

Finally, £ G L 2 (R + ) implies cu = 0. □ 

Recall the construction of the "forbidden" subspaces c C N , a ^ 1, given in 
[CK06b] . Let V = V* G £ 2 ([0, 1]; C WxiV ). For each a^l denote 

= [5 Q (^)] ± , where S a = 5 a (F) = f\(p*(p](t, A Q , V)dt = 5* > 

Jo 

and £ a = RanP a . Note that dimjFo, = N — dimS a = N — k a . The main result of 
[CK06b] is that one can modify each particular projector P a (keeping the spectrum 
and all other projectors fixed) in an arbitrary way such that T a D RanP a = {0}. It's 
quite natural that this restriction is equivalent to property (C) as shows 

Proposition A. 4 (Connection between subspaces T a and property (C)). 

Let (3^1 and (X a ;P a )^ L = (X a (V) ; P a (V))+^ for some V = V* G £ 2 ([0, 1]; C£ xjv ). 
Then, the collection (A a ; P a )a=i! where P a = P a for all a 7^ (3, satisfies (C) iff 

Tp n RanP^ = {0}, where Tp = [Sp(£p)] L . (A.4) 

Moreover, Tp = [Ran^ /3 (A ( g)]- L ; where £p is given by l(A.3\) . 

Proof. It follows from Proposition IA.3I (ii) that (C) holds true for the new collection 
(A Q ; Pa)t=i if and only if Pp£p(Xp)h ^ for all h G £% h 0. In other words, (C) is 
equivalent to 

Rane /3 (A /3 )nKerP /3 = {0}. (A.5) 
One has (see Lemmas 2.4 and 2.1 [CK06bJ for details) 

£p(Xp) = x(0, A^)Pj = Xp)P* = -<p*(l, XfiWl, Xp) X '(0, A^)P». 

Moreover, Ranx'(0, Xp)Pl = £p and 

Ran^g(A^) = Ran[0*(l, \p)ip'(l, \p)Pp] = RanSpPp = Sp(£p). 

Since dim Ker Pp = N - kp = N - dim Sp(£p), (TA5l) is equivalent to (TA4|) . □ 
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We finish our discussion by the consideration of the special case when only finite 
number of P a differ from the standard unperturbed coordinate projectors. 

Let A = {«!, a 2 , .., a m } be some finite set of exceptional indices. Assume 
that P a = P° coincides with some coordinate projector P° for all a ^ A 
(we admit multiple eigenvalues). Introduce the sets 

A] = {a £A : P a e] ± 0} 

(possible multiple eigenvalues belong to several A®). Assume that there 
exists C > such that the set {\ a ,a G A®} fl (— oo,7r 2 n 2 + C] consists 
of exactly n—m points for all j = 1, 2, .., N, if n is large enough. Let 

k ai + k a2 + . . + k am = Nm 

We give the simple description of all finite sequences (P a3 )T=n r ankP a = k a , such that 
the whole collection {(X a ; P a )}t=i satisfies (C): 

Proposition A. 5. Let (\ a ',Pa)t=i be as described above. Then (C) holds true iff 

( T Ti ... T m _i \ 
Ti T 2 ... T m 



T = 



T 



= T* > 0, 



2m-2 / 



where 



T k = Yl A*F(A a )P a F(A Q ) =T*, k = 0, 1, ., 2m-2, 



aeA 



F(A) = diag{/i(A),/ 2 (A), ..,/*(>)} and /,(A) = J[ 



eA° 



Remark A. 6. Since T)0in any case, the condition T > is equivalent to detT ^ 0. 

Proof. Indeed, let f(A) = (£i(A), 6(A), .., &v(A)) T be such that £(A) = 0(el Im ^l), 
£ G £ 2 (K+) and P a l(K) = for all a ^ 1. In particular, P°£(A Q ) = for all a G 
In order words, z^(z 2 ) G PW^_i 5 i] and £j(Aq,) = for all a G A®. Therefore, 

0(A) = Qj{\)fj{\), degQj^m-1, 
for some polynomials Qj. Let 

m—1 

g(A) = (g 1 (A),g 2 (A),..,g w (A)) T = ^A%, y p GC w and y = (y p )™^ G C 7 



iNm 



Then, 



m—1 



p=0 



m—1 



p,g=0 



p,<3=0 L a£A 



y 9 



= E [g(A a )]*p(A a )p a p(A a )g(A a ) = iaK)YP a aK). 

aeA aeA 

Hence, the NmxNm matrix T is degenerate iff there exists £ such that P Q £(A Q ) = 
for all a G A (recall that P^^Aq.) = holds true for all a ^ A by the construction). □ 
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B. Appendix. Three classical choices of additional spectral data in 

THE SCALAR CASE. 

In the scalar case, it is well known that the Dirichlet spectrum a(q) = {A n (g)}^^ 
determines only "one half of the potential q. Thus, in other to determine q uniquely, 
one needs either to assume that some partial information about q is known or to 
consider some additional spectral data besides o~(q). Note that there are two classical 
assumptions about the potential that make the knowledge of the spectrum sufficient: 
symmetry q(x) = q(l — x) (see, e.g., [PT87] ) or the knowledge of q(x) as a; 6 [0, |] 
(the Hochstadt-Lieberman theorem [HL78j . see also jGSOOj . jHo05] . [MP05] k Also, 
there are several classical choices of additional spectral data: 

(1) The second spectrum. This setup goes back to the original paper of Borg |Bo46j . 
The most natural choice is the spectrum {/U ra (g)} n= i of the mixed problem 

-y" + qy = Ay, y(o) = y'ii) = o. 

Note that {fi n (<l)}n=i U {^n(l)}n=i i s the Dirichlet spectrum of the symmetric 
potential q{2 — x) = q(x), x G [0, 1], defined on the doubled interval [0, 2]. 

(2) The normalizing constants (firstly appeared in Marchenko's paper [Mar50j) 

[ttn(g)]" 1 = / (f 2 (x,X n )dx = [yV] -1 ^, A n ) = - ^tttv^ = ~ res m W- 

JO J Xl U , A n ) A=A n 

(3) The norming constants introduced by Trubowitz and co-authors (see [PT87j ) 

-i y(-A) ' 

It is quite well known in the folklore that the characterization problems in the setups 
(l)-(3) are equivalent. Unfortunately, we do not know the good reference for this fact. 
So, the main purpose of this Appendix is to give the short proof of these equivalences 
(note that our arguments are quite similar to |Lev64] ) . For the simplicity, we assume 
that q G L 2 (0, 1), J q(x)dx = 0, i.e., {A n (g) — 7r 2 n 2 }^ G £ 2 (the similar arguments 
work well for other classes of potentials and corresponding classes of spectral data). 
Note that 

Hi < Ai < fi 2 < A2 < 1^2 < ■■■ an d /i n = Ti 2 {n — |) 2 + 0(1) as n — > 00. (B.l) 

Also, the Hadamard factorization implies 

+00 . . +00 . 

/(A)=<p(i,A)=n^Fr and g(A)=yUA)=n ^ 

- LJ - TT Z m Z - LJ - TT Z [m+i:) Z 

m=l m=l v 

Recall that we write a n = b n + £\ iff {n k \a n -b n \} n ~=i G £ 2 . 

Proposition B.l. Let X n = 7r 2 n 2 + £ 2 , ^B. 1\) hold and f(X), g(X) be given by (ELlf- 
Then, the following conditions are equivalent: 

(1) The asymptotics fi n = 7i 2 (n — |) 2 + £ 2 hold true. 

(2) The asymptotics a n = g(X n )f(X n ) = (27r 2 n 2 ) -1 (l + £\) hold true. 

(3) The asymptotics v n = log[(— l) n ^(A n )] = if hold true. 



i/ n (g) = log[(-l)V(l,An)] = log 



:-i) r 
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Proof. We start with the equivalence (2) <^> (3). Denote X n = n 2 X n — n 2 = 0(1) as 
n — ► oo. Then 



/(An) = 

Note that 



7T 2 72 2 



n 



Ajn A n 

7r 2 m 2 



( — 1)™ t-t A m — A n _ (— l) n 
2ir 2 n 2 J-J- n 2 (m 2 — n 2 ) 2-K 2 n 2 



n 



i + 



Am, A n 

m 2 — n 2 



log J] 

m^n 



1 + 



Am A n 

m 2 —n 2 



Am A n 

m 2 — n 2 



+ 



m^n L 

A m — A n ^ ^ 

m^ra 



(m 2 —n 2 ) 2 



m?—n 2 



n- 



E 

m^n 



m? — n 2 



+ 



rr 



Then, it immediately follows from (A n )^^ G £ 2 and simple properties of the discrete 
Hilbert transform (see Lemma E2 (ii) below) that /(A„) = (-l)"(27r 2 n 2 )- 1 (l 
Thus, (2) <^ (3). The proof of the equivalence (1) <^ (3) is similar. Indeed, 



+oo 



where fi r . 



H"m A n 

- 7r 2 (m — \) 2 



+oo 



A'm A n 



+oo 



- 7T 2 ((m — |) 2 — n 2 ) 



m=l L 



1 + 



fJ"m A n 



(m — |) 2 — n 2 



IX 2 /i r 



m + tj) 2 = 0(1) as m — > oo. As above, 



+oo 



log[(-l)^(A n )] = 



/^m A n 

v m — i) 2 — n 2 

m=l v ^' 



v 7 m=l 



(m — |) 2 — n 2 



+ — 



and the equivalence (1) (3) follows by Lemma IB. 21 (i). 
Lemma B.2. (i) The linear operator (a m )m=i {b n )t=i> where 

1 +oo 1 + 00 r- 

\wn ^ ffi _( m — 1)2 tt 



□ 



27m 



m=l 



m=l 



n— m + i n — (l—m)-\-\ 



oo 

m)m=l 



is an isometry in £ 2 . 
(ii) The linear operator (a 

is bounded in £ 2 . 



(K)t=i> where 

+00 



2n ^— ' n 2 — m 2 

m=l 



E 

m=l 



a.„ 



— hi, | 

n — m n—[—m) 



Proof. Both results easily follows by the Fourier transform and the identities (in L 2 (T)) 



E 



Til 



fc=— oo ^ 

where C = ^ 1, <j> E (0,2vr). 



7rie~^~ 



and 



E 



-2 



7T 



□ 



Remark B.3. TTie similar technique can be applied for other inverse problems in order 
to derive the characterization of some additional spectral parameters (e.g., similar to 
a n (q) ) from the characterization of other parameters (e.g., similar to u n (q) ). In general, 
these characterizations may differ from each other substantially, see [CK07J . 
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